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Abstract 

The probability that interferometric detectors such as LIGO and VIRGO will suc- 
cessfully detect inspiraling compact binaries depends in part on our knowledge of the 
expected gravitational wave forms. 
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ff^ , The best approximations to the true wave forms available today are the post- 

^^ ' Newtonian (PN) templates. In this paper we argue that these 2PN templates are 

I . accurate enough for a successful search for compact binaries with the advanced LIGO 

t^^ I interferometer. Results are presented for the 40-meter Caltech prototype as well as for 

^ ' the inital and advance LIGO detectors. 
O 

^: 1 INTRODUCTION 

Currently (November 1997) there are five large scale interferometric gravitational-wave de- 
rS ' tectors under construction, namely the two American LIGO detectors, the French Italian 

c^ ■ VIRGO, the German-British GEO600 and finally the Japanese TAMA. The immediate goal 

of this effort is the direct detection of gravitational waves. In a later stage this international 
network of detectors will act as a gravitational wave observatory with a variety of scientific 
projects (see reference |l[] for more details). The most promising sources for detection are 
inspiraling binaries of black holes and/or neutron stars. 

Because since gravitational waves are so weak, any signal will likely be obscured by 
detector noise. The method of choice to search for a known signal in a noisy data stream is 
matched filtering 0: The data is compared to a set of template wave forms for a "best fit" (a 
more precise definition is given below). Matched filtering is an optimal method in the sense 
that no other (linear) filtering method is more effective. The drawback is that a set of very 
accurate templates is needed to track the signal with the necessary accuracy. For example 
a typical inspiraling-binary signal might approximately 15000 wave cycles in the advanced 
LIGO frequency band. This means that we need a template that traces all of the 15000 
cycles. By missing only one cycle the detection probability falls drastically. Unfortunately, 
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in general relativity there is no way to (analytically) calculate the exact wave-form and we 
have to resort to some kind of approximation. The most promising approximative templates 
are obtained by a post-Newtonian expansion of the gravitational waves in the orbital velocity 
of the binary companions. 

In this contribution we assess the quality of the PN templates by comparing them to 
signals obtained from black hole perturbation theory. This is known to provide very accurate 
wave forms when the binary mass ratio is small. 

2 MATCHED FILTERING 

Let us write the detector output as 

s{t) = h{t)+n{t), 

where n{t) denotes the detector noise and h{t) a possibly absent signal. We define the 
squared signal to noise ratio (SNR) by 
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Here the inner product ( , ) is defined by 
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where a hat "denotes a Fourier transform and an asterisk * complex conjugation. S'„ is the 
(one sided) spectral density of the noise. It is large where the noise is large. Consequently, 
regions of low noise get a higher weight in the integral (^), the "match" between the template 
t and the detector output h. If in an actual measurement p exceeds a certain threshold then 
detection is claimed. 

Wieners theorem f^ now states that p is maximal if t = /i in equation (|l]). This choice 
of template is called the Wiener optimal filter. Conversely, ii t ^ h then we get a smaller p 
and thus might miss signals that otherwise would have been detected. 

Since we want to make statements about the templates for "typical" noise, we usually 
talk about the the expectation value of the SNR. Under the assumption of Gaussian noise it 
can be calculate to be [§ 
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The quantity (p )best = y{h, h) is the SNR one would obtain if one knew the exact waveform 
h. The function A is called the ambiguity function and it depends on the template used. 
Maximizing A over all templates gives the fitting factor FF . This is the reduction of the 
SNR caused by inaccurate templates: (p^)max = FF{p^)\,cs\.- The fitting factor therefore 
measures the quality of the templates. 



The best templates available today are calculated using a post-Newtonian expansion 
of Einstein's equations. For inspiraling binaries, the wave forms are expressed as a power 
series in v/c where v is the orbital velocity. (Recently Fade approximates have been used to 
improve the convergence of the post-Newtonian expansion P] ) . 

To simplify matters, one often works with the so-called restricted post-Newtonian ap- 
proximation given by 

t(/) oc /-Ie^*(-^.'?.*c,</.c)^ (3) 

where the amphtude is kept at the lowest post-Newtonian order but the phase \1/(A^, r^, tc, 0c) 
includes higher FN corrections. Here we keep \1/ accurate to 2 FN order. This makes sense, 
because we expect a slight mismatch in the phases to lead to a rapidly oscillating integrand 
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in (^ and thus to a low SNR. Here the chirp mass M. = {mim2)'^ /{mi + 1712)'^, the mass 
ratio 1] = {17111712) / [mi + 1712)^, the arrival time tc and the initial phase 0c parameterize the 
templates. 

The test signal used to calculate the fitting factor is obtained from black hole perturbation 
theory |^, integrating the Teukolsky equation 0. The signal is given as an infinite sum over 
modes and takes the form 

M/)ocEE^'"^Me'^^^^5K^,^) (4) 

1=2 m=l 

where v = (27rM//'m)^/^ (M = 1711 + 1712)- The functions Aim{v) and 'ip{v) must be calculated 
numerically and Si is a known function. The I = m = 2 mode dominates the signal and we 
have verified numerically that modes with / > 3 can be neglected in the sum (^). 

To calculate the fitting factor one sticks expressions (^ and (^ into the definition of the 
ambiguity function and maximizes over the parameters A4, i], tc and 0c for a "fixed" signal 



3 RESULTS AND INTERPRETATION 



advanced LIGO initial LIGO 40-meter 

System FF-Newton FF- 2-FN FF- 2-FN FF- 2-FN 



1.4- 


- 1.4M0 


79.2% 


93.0% 


95.8% 


94.0% 


0.5- 


- 5.OM0 


51.6% 


95.2% 


89.7% 


89.0% 


1.4- 


-IOM0 


55.8% 


91.4% 


88.4% 


< 85.0% 


10- 


-IOM0 


70.1% 


91.7% 


86.3% 




4- 


■3OM0 


61.3% 


86.6% 


67.9% 





Table 1: The fitting factors for different noise curves and binary systems are listed. A fitting 
factor larger than 90% is considered acceptable. Note that only the 1.4 — I.4M0 system, for 
which FF > 90%, is within the targeted range of the 40-meter detector. 

We have calculated the fitting factors for a variety of different binaries and detector noise 
curves (for the initial and advanced LIGO detectors as well as the 40m Caltech prototype). 



In the former case, in addition to the 2 PN templates we have also tested the Newtonian (= 
lowest PN) templates. The results are summarized in Table 1. 

We find that for the mass range of interest to a particular detector, the fitting factors for 
the 2PN templates are all above 90%. This makes them suitable for detection of inspiraling 
binary signal |^. For small mass ratios the fitting factors drop considerably refiecting the 
fact that the systems become more relativistic for fixed chirp mass Aic- One would expect 
PN templates to perform less well in this limit. 

An additional point should be made. Formally, black hole perturbation theory assumes 
rj = 0. Thus it is not clear if our analysis holds for nearly-equal-mass systems. However, by 
setting ?7 = in the templates (0) one effectively obtains a PN expansion of the perturbation- 
theory signal. When we do so we do not find any significant changes in the fitting factor. 
This leads us to believe that our results should hold, at least qualitatively, even for nonzero 
mass ratios. 

In conclusion, we may say that the 2 PN templates are suitable for the detection of binary 
inspiral signals. 
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